To characterize the molecular environment of classical 
T Tauri stars in Lupus observed with the Spitzer Space 
Telescope and to search for gas-rich disks toward these sources. 
Submillimeter observations of '-CO, '^CO and C'^O in the 
7=3-2 and 2-1 lines were obtained with the James Clerk 
Maxwell Telescope toward a sample of 21 T Tauri stars with 
disks in the Lupus molecular clouds. Pointings at the sources 
and at selected off-positions are presented in order to disentan- 
gle material associated with disks from ambient cloud mate- 
rial. One source, IM Lup (Sz 82), has been found with double- 
peaked '^CO and '^CO profiles characteristic of a large rotat- 
ing gas disk. The inclination of the disk is found to be ~ 20°, 
with an outer radius of 400-700 AU. For most other sources, 
' including GQ Lup with its substellar companion, the single- 
' dish '^CO lines are dominated by extended cloud emission 
with a complex velocity structure. No evidence for molecular 
I outflows is found. Compact C'^O emission due to dense cir- 
P5 cumstellar material has been detected toward only two sources. 
^ , Future searches for gas-rich disks in Lupus should either use 
• interferometers or perform very deep single-dish integrations 
' in dense gas tracers to separate the cloud and disk emission. 

Key words. Star formation - low-mass stars - Lupus - T Tauri 
J> star - circumstellar disk - IM Lup - GQ Lup 
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Abstract. 



1. Introduction 

A general characteristic of pre-main sequence (PMS) stars 
is the presence of circumstellar disks, which have been ob- 
served for a wide range of evolutionary stages associated with 
low-mass star formation (e.g., Lada 1987, Andre et al. 1993, 
Greaves 2005). These disks contain all the necessary ingredi- 
ents for forming complex planetary systems similar to our own 
Solar system, now also seen around other stars (e.g., Ruden 
1999, Marcy & Butler 1995). The gas-rich phase of disks is 
crucial in the evolution, since gas is needed to build Jovian- 
type planets. Also, the presence of gas affects the dynamics 
and evolution of the dust in the system and the migration of 
any planets embedded in the disk (e.g., Gorti & Hollenbach 
2004). With time, the disk will loose its gas through processes 
such as photoevaporation, stellar winds, formation of gas-rich 
planets and accretion onto the star (e.g., Hollenbach et al. 2000, 
Alexander et al. 2005). However, the timescale for gas dissipa- 
tion in disks is still poorly constrained observationally and it is 
even unclear whether gas and dust are lost simultaneously from 
disks when they evolve from the massive optically thick to the 
tenuous optically thin 'debris' phase. 

To constrain these timescales, surveys of gas in disks in a 
variety of environments are needed with sensitivities sufficient 
to detect a Jupiter mass of gas (~10"^ M©). Surveys to date 
have focussed on CO milUmeter observations, mostly of disks 
in the Taurus molecular cloud (e.g. Thi et al. 2001). Although 
CO is known to be a poor tracer of the gas mass due to the 
combined effects of photodissociation in the upper layers and 
freeze-out in the midplane (e.g., van ZadeUioff et al. 2003), it 
is much simpler to observe with current instrumentation than 
other tracers including H2 (Thi et al. 2001) or atomic fine- 
structure fines (Kamp et al. 2003). Interferometer surveys have 
detected gas-rich disks around a large fraction (>60%) of clas- 
sical T Tauri stars in Taurus- Auriga with ages up to a few Myr 
and gas masses of at least a few times 10"^ Mq, depending 
on assumptions about the CO abundance (Koerner & Sargent 
1995, Dutrey et al. 1996; 2003). In contrast, the detection rate is 
only 10% for the more evolved weak-line T Tauri stars (Duvert 
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et al. 2000). Similar statistics are found in a single-dish submil- 
limeter survey of CO in a spatially more distributed sample of 
Herbig Ae and Vega excess stars (Dent et al. 2005) and in CO 
infrared surveys (Najita et al. 2003, Blake & Boogert 2004). 

Surprisingly fittle is known about the presence of gas-rich 
disks around T Tauri stars in other nearby star-forming re- 
gions such as Corona Australis, Lupus and Chamaeleon. This 
is largely due to the absence of interferometers in the Southern 
sky, and even Ophiuchus has been poorly sampled with current 
facilities. We present here an initial search for gas-rich disks in 
the Lupus molecular clouds, using a single dish telescope. 

The Lupus clouds, located at around galactic coordinates 
335° < ; < 341° and T < b < 17° (Hughes et al. 1994), 
are among the closest star-forming regions to the Sun at ~ 150 
pc (see reviews by Krautter & Reipurth 1991, Comeron 2006). 
Mid M-type stars dominate the stellar population of Lupus, 
which contains no confirmed O or B-type stars. In total 69 
Pre-Main Sequence (PMS) objects were found by Schwartz 
(1977), distributed over four clouds now known as Lupus 1^ 
(see also Hughes et al. 1994). In addition, some 130 new weak- 
line T Tauri stars were detected in ROSAT images (Krautter 
et al. 1994; 1997). The most massive stars in the complex are 
two A-type stars, HR 5999 and HR 6000, present in the rich- 
est sub-grouping in stars and complexity. Lupus 3. A good 
overview of the clouds is given by the extinction map made 
by Cambresy (1999), based on star counts. Recently a higher 
resolution (~ 30") map has been pubUshed by Teixeira et al. 
(2005). 

At milfimeter wavelengths the Lupus clouds have been 
studied both in molecular emission lines and in the continuum. 
i^CO (Tachihara et al. 2001), '^CO (Tachihara et al. 1996) and 
C'^O (Hara et al. 1999) 7=1-0 large-scale maps were made at 
2^6 resolution with the NANTEN telescope. Higher resolution 
observations were done with the Swedish-ESO Submillimeter 
Telescope (SEST) over more limited regions with ~45" reso- 
lution by Gahm et al. (1993) (mostly '^CO 1-0 in Lupus 2), 
Rizzo et al. (1998) ('^CO 1-0 in Lupus 1 and 4 filaments) and 
Vilas-Boas et al. (2000) ('^CO and C'^^O 1-0 in locaUzed dark 
cores in Lupus 1-4). Niirnberger et al. (1997) performed a con- 
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tinuum survey of 32 T Tauri stars in Lupus at 1.3 mm with the 
SEST, probing the cold circumstellar dust. 

The Lupus star-forming region appears to differ in several 
aspects from other low-mass star-forming clouds. In contrast 
with Ophiuchus and Taurus, which contain a significant frac- 
tion of deeply embedded Class and Class I objects, Lupus 
has at most a few embedded objects (Krautter & Reipurth 1 99 1 , 
Tachihara et al. 1996, Comeron 2006). This could be an indi- 
cation that Lupus is more evolved than other star-forming re- 
gions. The peak in the age distribution is 3 x 10^ yr for stars in 
Lupus 3 and 4, calculated from the PMS evolutionary tracks of 
D' Antona & Mazzitelli (1994) for a distance of 150 pc (Hughes 
et al. 1994). Under the same assumptions. Lupus 1 and 2 PMS 
stars are younger with a peak of 1 x 10^ yr. These age esti- 
mates are comparable to those for other star-forming regions 
such as Chamaeleon (Hartigan 1993), but are about double that 
of the Taurus cloud (Simon et al. 1993). Thus, a study of disks 
in Lupus is interesting because this is precisely the age range 
in which significant evolution of the gas disk is expected. 

A second characteristic of the Lupus clouds is the distribu- 
tion of spectral types. Compared to the Taurus cloud, Lupus is 
dominated by lower mass stars. The peak of the distribution in 
spectral type is MO, with very few stars having spectral types 
higher than K7. In contrast, the peak of the spectral type distri- 
bution in Taurus is at K7. 

In this work, single-dish submillimeter observations of CO 
toward a sample of 21 classical T Tauri stars with disks are pre- 
sented obtained with a beam of 14", much smaller than that of 
previous data. The sources are a subset of the PMS stars found 
by Schwartz (1977), Krautter & Reipurth (1991) and Hughes 
et al. (1994), and they have been observed with the Spitzer 
Space Telescope in the context of the 'From Molecular Cores to 
Planet Forming Disks' (c2d) Legacy survey (Evans et al. 2003). 
The aim of the CO observations is to characterize the molecu- 
lar environment of the T Tauri stars at high angular resolution 
and to search for gas-rich disks suitable for future follow-up 
interferometer observations. 

2. Observations 

2.1. CO observations 

The submillimeter CO observations were carried out in two 
runs in April and July 2005 at the James Clerk Maxwell 
Telescope (JCMT)'. '^co 7=3-2 and C'^^O 7=3-2 and 2-1 
transitions (depending on weather conditions) were observed 
with the dual polarization B3 and single polarization A3 re- 
ceivers. The beam sizes are 14" and 20", respectively. The 
Digital Autocorrelator Spectrometer was used as the back-end 
with a bandwidth of 125 MHz, giving a spectral resolution of 
98 kHz for the A-receiver corresponding to 0.10 km s~' at 230 
GHz, and 196 kHz for the dual-polarization B-receiver corre- 
sponding to 0.22 km s~^ at 345 GHz. The integration times 

' The JCMT is operated at the summit of Mauna Kea, Hawaii 
by the Joint Astronomy Center on behalf of the United Kingdom 
Particle Physics and Astronomy Research Council, the Netherlands 
Organization for Scientific Research, and the National Research 
Council of Canada 



were such that the RMS is around 70-100 mK in a 0.2 km s~' 
velocity bin for the B-receiver and 50 mK for the A-receiver, 
depending on source and observing conditions, on a T^h scale. 
Pointing was checked regularly for these southern sources and 
found to vary within 2". Beam efficiencies were taken to be 
0.65 for RxA3 and 0.63 for R\B3^. 

Frequency switching with a switch of 16.4 MHz was 
adopted for ^^CO because no emission-free position could 
readily be found for many of the sources in our sample. If a sig- 
nal was detected on source, additional observations were car- 
ried out at 30" offsets to the east and south of the source to 
characterize its environment and constrain the on-source con- 
tribution. If indications were found of a gas-rich disk, addi- 
tional follow-up observations were done on-source in ^^CO 3- 
2 using a beam switch of 180" in AzEl. C'^O 3-2 or 2-1 ob- 
servations were performed for nearly all sources, also using a 
beam switch of 180". Table 1 summarizes the observed lines 
per source. The data were reduced with the SPECX and CLASS 
reduction packages. 

Our observing strategy is similar to that for the single-dish 
observations of gas-rich disks in Taurus by Thi et al. (2001). 
Because the surface layers of disks are warm, the '^CO 3-2 line 
is well suited to observe these regions. Also, the smaller beam 
size at higher frequencies is more favorable for observing disks, 
which have typical sizes of a few hundred AU corresponding 
to a few arcsec at the distance of Lupus. Typical ^^CO 3-2 an- 
tenna temperatures of disks in Taurus are 0.5-1 K, so that an 
RMS of <0. 1 K needs to be reached. To distinguish gas asso- 
ciated with disks from that of surrounding clouds, observations 
on-source and at nearby off-source positions need to be taken. 
Gas-rich disks should show a characteristic double-peaked line 
profile and a stronger signal on-source than at the off- source 
positions, while cloud emission should produce a comparable 
signal at all positions. Finally, disks in Taurus have not been de- 
tected in C^^O with single-dish telescopes. On the other hand, 
such beam-switched C'**0 observations have been found to be 
a good probe of any dense circumstellar gas or remnant enve- 
lope associated with the source. For example. Class I sources 
show typical C'^O intensities of 0.5-2 K (Hogerheijde et al. 
1998, J0rgensen et al. 2002). 

2.2. Sample 

Our sample of classical T Tauri stars in Lupus was selected 
from targets belonging to the InfraRed Spectrometer (IRS) 
sample of the c2d Legacy program on Spitzer (See Table 1). 
They are a subset of the T Tauri stars studied by Schwartz 
(1977), Krautter & Reipurth (1991) and Hughes et al. (1994), 
selected to be brighter than 200 mJy at 12 fim and excluding 
sources contained in the Spitzer-IRS guaranteed time programs. 
Infrared photometry and spectroscopy was performed on all 
sources giving estimates of their spectral types (see Table 1). 
In total 21 PMS stars were observed, most of them located in 
Lupus 3. All sources show a clear 10 jum emission feature due 
to warm siUcates, confirming the presence of a circumstellar 

^ Observed by the JCMT staff over many observations; see the 
JCMT website http://www.jach.hawaii.edu/JCMT/ 
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disk (Kessler-Silacci et al. 2006). Some of the selected stars 
are known binaries. Where available, results on the 1.3 mm 
continuum emission are reported (Ntirnberger et al. 1997). This 
emission is presumably dominated by the cold dust in the disk. 

The distances to individual Lupus clouds have consid- 
erable uncertainty. Hughes et al. (1993) found 140±20 pc 
from spectroscopic parallax of field stars around the clouds. 
However, observations of reddening and interstellar absorption 
lines toward individual stars by e.g., Knude & H0g (1998) and 
Crawford (2000) find distances ranging from 100 to 200 pc. A 
good discussion can be found in Comeron (2006). The depth 
of the Lupus clouds may be non-neglible and can be as large 
as 50 pc. Also, it has been suggested that not all clouds of 
Lupus are related, especially Lupus 2 could be located as far 
as 360 pc (Knude & Nielsen 2001). However this seems un- 
likely given the narrow ranges in mean velocity of the clouds 
(see e.g., Tachihara et al. 1996, Vilas-Boas et al. 2000 and this 
work §3). The distance used here for all sources is that most 
commonly adopted, 150 pc. 

3. Results 

3.1. i^CO 

The observed ^^CO 3-2 Une intensities are presented in Table 2, 
together with the peak intensities. Typical examples of spectra, 
illustrating the range of observed profiles, both on-source as 
well as at off-source positions, are shown in Fig. 1. The re- 
mainder of the spectra can be found in the on-line material. 

Only a single source is found in our sample with a clear 
signal of gas associated with the source itself, namely IM Lup. 
A classic double-peaked disk profile is found toward the source 
position, which is largely absent at the 30" east and south po- 
sitions offset from the source (see Fig. 1). The profile, with a 
maximum peak of 2.5 K, is significantly broader than that as- 
sociated with quiescent gas. This source is further discussed in 
§4.1. 

None of the other sources show indications of gas-rich 
disks surrounding the central T Tauri star. Three sources, RY 
Lup, EX Lup and Sz 111, appear to be located outside any 
cloud material and have no detectable disk emission at our sen- 
sitivity limit. Using the formula in Thi et al. (2001) with the ob- 
served limits on ^^CO emission, upper limits of the gas masses 
of 1 X 10"^ Mq are obtained for these disks with a CO/H2 ratio 
of 10""*. For a more realistic overall CO/H2 disk abundance of 
10"*, the limits would increase to 1 x 10""* (see also §4). 

Most sources have a strong '^CO line on source, but these 
signals are also detected at the off-positions with comparable 
strengths (see Fig. 1). These lines must originate in the large 
scale molecular cloud material in the line of sight of the source. 
The complex profiles seen for many sources indicate a com- 
plex structure of different molecular clouds moving at different 
velocities, either in the fore- or background. For example, the 
spectra toward Sz 1 18 (see Fig. 1) show three peaks with differ- 
ent intensities at all positions with widths of ~0.8 km s~^ and 
velocities that differ by 1-1.5 km s"'. C^^O is not detected at 
all for this source. Comparison with the large scale maps shows 
that observations of ^^CO 3-2 with peak temperatures of 2-20 



Kelvin are consistent with the lower resolution '^CO surveys 
(Gahm et al. 1993, Tachihara et al. 1996). 

The i^co fine widths of -0.8 km s ' are comparable to 
those found in Vilas-Boas et al. (2000) and Tachihara et al. 
(2001), taking into account multiple clouds at different veloci- 
ties. No evidence for line wings at >2 km s~^ from Une center is 
found at the level Tub < 0.15 K (2<t), providing limits on any 
small-scale molecular outflows associated with these T Tauri 
stars. For RU Lup, this is consistent with Gahm et al. (1993), 
who reported no outflows from CO 1-0 data in a larger beam. 

Although no direct correlation can be inferred due to the 
confusion of cloud emission around most sources, it is interest- 
ing to note that the largest 1.3 mm flux found by Niirnberger 
et al. (1997) is associated with IM Lup. Thus it also contains 
the most surrounding dust, in addition to being the only source 
with CO disk emission. 

3.2. C^^Oand^^CO 

The results for the C'^0 and '^CO line intensities are included 
in Table 2. '^CO 3-2 was observed for only two objects: IM 
Lup and GQ Lup. A more detailed description for these two 
objects can be found in §4. Figure 2 shows some examples of 
observed C'^O lines. 

Only two sources in our sample show a C^^O fine above the 
noise level: HT Lup and Sz 73. HT Lup is the only source to- 
ward which both C'**0 2-1 and 3-2 lines have been detected 
in our survey. Because the data were taken in 180" beam- 
switched mode, the bulk of the emission is likely directly as- 
sociated with the source. This could point towards a small cir- 
cumstellar envelope. The 2-1/3-2 line ratio of > 1 indicates 
either very cold gas (< 15 K) or densities less than 10^ cm"^ 
for r s; 20 K. The inferred C'**0 column density for these con- 
ditions is ~ 10^^ cm"^, corresponding to a CO column density 
of 5 X 10" cm-2 for i«O/i»O=540 (Wilson & Rood 1994) and 
an H2 column density of 6 x 10^^ cm"^ for a dark cloud C'^O 
abundance of 1.7 x 10"^ (Frerking et al. 1982). Thus, the ex- 
tinction directly associated with the source is Ay « 6 mag using 
Nh/Av = 1.8 X 10^1 cm-2 mag"' (Rachford et al. 2002). 

As mentioned before, our observations are not sensitive to 
extended C'^O emission because of beam-switching by 3'. The 
lack of any features at the 2cr level of <0.1 K in C'**0 2-1 and 
<0.3 K in 3-2 implies that either any large scale cloud emission 
must be very smooth on 3' scale or that there is no C'^O at 
this limit. Vilas-Boas et al. (2000) and Hara et al. (1999) find 
C'^O 1-0 antenna temperatures of 0.2 to 1.5 K for the densest 
cores. For typical dark cloud conditions, the 2-1 line should be 
comparable in strength whereas the 3-2 line should be factors 
of 2-4 weaker, so that these lines would have been detectable at 
the positions of strongest 1-0 emission. However, most of our 
sources are located at the edge or outside the C'^O contours 
of Hara et al. (1999) suggesting that C'^O is indeed largely 
absent. The lack of C^^O 2-1 emission at tiie level of 0.08 K 
(2cr) indicates a C'^O column density of less than 5xlO'^ cm"-^, 
or H2 less than 3 X lO^'' cm"^ using the dark cloud abundance. 
In this diffuse cloud regime, however, the C^^O abundance may 
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Fig. 1. '^CO 3-2 lines of selected sources in Lupus, with off-positions 30" south and east of source. 



be significantly lower than the standard dark cloud value (van 4. Individual Objects 
Dishoeck & Black 1988), increasing the H2 column density. 

4.1. IMLup 



IM Lup is the only T Tauri star in our sample with a clear de- 
tection of ^^CO emission associated with the source. A double- 
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Table 1. Summary of the properties of T Tauri stars and observed lines in Lupus 



Source 


RA (J2000) 


Dec (J2000) 


Spectr.type" 


Cloud 


^1 .3mm (mJy)'' 




Offsets 


c'«o 


Add. obs 


HT Lup" 


15 45 12.87 


-34 17 30.6 


K2 


Lupus 1 


135 


3-2 


E S 


3-2/2-1 




GWLup 


15 46 44.68 


-34 30 35.4 


M2-M4 


Lupus 1 


106 


3-2 


ES 


3-2/2-1 




HMLup 


15 47 50.63 


-35 28 35.4 


M3 


Lupus 1 


<45 


3-2 


ES 


3-2/2-1 




Sz 73 


15 47 56.98 


-35 14 35.1 


MO 


Lupus 1 


26 


3-2 


ES 


2-1 




GQ Lup 


15 49 12.12 


-35 39 05.0 


K7-M0 


Lupus 1 


38 


3-2 


E S 


2-1 


"CO 3-2 


IMLup 


15 56 09.17 


-37 56 06.4 


MO 


Lupus 2 


260 


3-2 


E S 


2-1 


'^CO 3-2 


RULup 


15 56 42.31 


- 37 49 15.5 


K7-M0 


Lupus 2 


197 


3-2 


ES 


3-2 




Sz 84 


15 58 02.52 


-37 36 02.8 


M5.5 


Lupus 2 


<36 


3-2 


ES 


2-1 




RYLup 


15 59 28.39 


-40 21 51.2 


K4 


Lupus 3 




3-2 




3-2 




EX Lup 


16 03 05.52 


-40 18 24.9 


MO 


Lupus 3 




3-2 




3-2 




HO Lup^ 


16 07 00.61 


-39 02 19.4 


Ml 


Lupus 3 


<24 


3-2 


ES 


3-2 




Sz 96 


16 08 12.64 


-39 08 33.3 


M1.5 


Lupus 3 


<45 


3-2 


E S 


3-2 




Krautter's star 


16 08 29.70 


-39 03 11.2 


KO 


Lupus 3 


<30 


3-2 


ES 


3-2 




Sz 107 


16 08 41.78 


-39 01 36.4 


M5.5 


Lupus 3 




3-2 


ES 






Sz 109 


16 08 48.17 


-39 04 19.0 


M5.5 


Lupus 3 


_ 


3-2 


ES 


3-2 




Sz 110 


16 08 51.56 


-39 03 17.5 


M2 


Lupus 3 




3-2 


ES 


3-2 




Szlll 


16 08 54.74 


-39 37 17.5 


M1.5 


Lupus 3 








3-2 




V908 Sco 


16 09 01.84 


-39 05 12.1 


M4 


Lupus 3 




3-2 


ES 


3-2 




Sz 117 


16 09 44.34 


-39 13 30.4 


M2 


Lupus 3 




3-2 


ES 


3-2 




Sz 118 


16 09 48.68 


-39 11 17.2 


K6 


Lupus 3 




3-2 


ES 


3-2 




Sz 123 


16 10 51.58 


-38 53 13.7 


Ml 


Lupus 3 




3-2 


ES 


3-2 





" Hughes et al. (1994) 

* Niimberger et al. (1997), where available 
^ Known binaries (Niimberger et al. 1997) 
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Fig. 2. C'^O lines toward selected sources in Lupus. 

peaked profile with peak separations of about 2 km s~^ is seen 
in both '^CO and '^CO, consistent with a rotating gas disk (see 
Fig. 3). In addition, the lack of C'^O emission over 0.09 K in 
a 0.25 km s~^ bin implies the absence of a significant contri- 
bution of a dense circumstellar envelope. The existence of a 
large disk around IM Lup has recently been shown in scattered 
light images at visible and near-infrared wavelengths with the 
Hubble Space Telescope (Padgett et al, priv. comm. Schneider 
et al, priv. conmi.). 

The '^CO peak intensity for IM Lup is comparable to that 
of the brighter sources with gas-rich disks in Taurus observed 
by Thi et al. (2001). The much weaker ambient cloud material 
seen at the south ofl'-position (see Fig. 1) modestly afl'ects the 




-5 5 ^ 10 15 

Vlsr [km s ] 

Fig. 3. '^CO 3-2, "CO 3-2 and C"*0 2-1 spectra observed toward IM 
Lup (Sz 82) shifted by +4, +2 and K, respectively. 



disk profile, however, especially the red peak at Vlsr = 5 km 
s~^ The blue peak at Vlsr = 3 km s~^ is not contaminated 
and is taken as the actual disk signal for both peaks. The '^CO 
peaks of 0.2 K are comparable to those for Taurus disks and do 
not seem to be aff'ected by ambient cloud material 

The observed '^CO/'^CO integrated intensity ratio of 6.0 
± 1.5 indicates '^CO and '^CO optical depths of 2.3 and 0.15, 
respectively, assuming that the excitation temperatures are the 
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Source 


j'TuBdV 




Peak Tub 






'"CO 3-2 "CO 3-2 


'^CO 


U U J—Z 


L. U Z-i 


"CO 3-2 




[Kkms '] [Kkms"'] 


[K] 


LJ^l 


IJ^J 


[K] 


HTLup 


26.5 


15 


1 ft 


1 c 




GWLup 


9.1 


6.7 


< U. to 


< u.us 




HMLup 


6.4 


5.2 


^ O 1 < 


^ n no 




Sz 73 


15.8 


9.0 




1 n 
l.U 




GQLup 


13.4 <0.14'^ 


5.8 




A^ 


<0.25'^ 


IMLup 


4.7 0.52 


2.4 




^A AA 


0.35 


RULup 






<0.27 






Sz84 


0.6 


0.8 




^A AA 




RYLup 


* 


<0.25 


<U.jl 






EXLup 


b 


< 0.16 


<0.4U 






HOLup 


1.3 


0.8 


^ A 1 '7 
< U.l / 






Sz 96 


25.1 


12 


^ A 1 
< O.lO 






Krautter's star 


17.5 


11.7 


<U.ZJ 






Sz 107 


18.3 


12 








Sz 109 


34.8 


26 


<0.15 






Sz 110 


26.1 


21 


<0.25 






Sz 111 


* 


<0.2 


<0.30 






V908 Sco 


22.6 


13.7 


<0.26 






Sz 117 


11.4 


5.5 


<0.22 






Sz 118 


14.3 


7.2 


<0.27 






Sz 123 


2.4 


1.6 


<0.20 






" If no line was detected, a 2cr™j(see § 2.1) limit in 6v = 0.25 km s ' is j 


given. 





* No detection. See peak value for RMS 
' see § 4.2 



same for both isotopes and constant throughout the disk. The 
uncertainties are dominated by calibration errors. If the peak 
temperatures of the blue wing are used (1.5 for ^^CO and 0.4 K 
for '^CO respectively), optical depths of 3 and 0.3 are found. 

Thi et al. (2001) give several methods to derive disk masses 
from different observations. Assuming that most of the '^CO 
3-2 emission is optically thin and originates in the cold gas of 
the outer disk with a single excitation temperature, the total gas 
mass is given by 



M„as = 3x10" 



'r,, +0.89 



'-COI/l'^COjWHz/'^CO 



60 



104 



(1) 



-16.02/r,, 



:)7 



T„hdV Mo 



,1 -e--" lyiOOpcj 

The H2/CO and '^CO/'^CO ratios have been set to 10"* and 
60, following standard conversion factors. This results in a de- 
rived gas mass of 2 x 10""* M© for IM Lup atd = 150 pc, using 
Tex - 30 K. For comparison, the disk mass found from the 1.3 
mm continuum flux (Ntirnberger et al. 1997) is 6 X 10"^ Mq as- 
suming a gas/dust ratio of 100, Tdust = 30 K and a dust opacity 
coefficient k of 0.01 cm^g~'. The factor of 300 between these 
two numbers is comparable to that found for other disks (e.g., 
Dutrey et al. 1996, Mannings & Sargent 2000, Thi et al. 2001) 
and is commonly explained by the fact that CO is not a good 
tracer of the gas mass due to the combined effects of photodis- 
sociation by steUar and interstellar UV radiation and freeze-out 
of CO onto the grains in cold regions of the disk (van Zadelhoff 
et al. 2001, Aikawa et al. 2002). 



The '^CO profile can be used to obtain estimates of two 
other important disk parameters, the disk size and inclination. 
Assuming optically thick ^^CO 3-2 emission and foUowing the 
method described in Dutrey et al. (1997): 



Jmbdv 



■ COSi 



(2) 



where and /?out are the inner and outer radii, 5v the turbulent 
velocity and p a geometrical factor. is the telescope beam 
on the sky and i the inclination, with degrees being face-on 
with respect to the Une of sight. For IM Lup, - AU, a 
mean disk excitation temperature Tex of 25, 40 and 50 K and 
{p6v) of 0.3 km s"' are chosen. 

The disk sizes calculated from Eq. (2) are given in 
/?out Vcosi. Scattered light images indicate that the inchnation 
is of order 10° to 30°, so that cos 1 is a reasonable approxi- 
mation. The disk sizes found with Tex = 25 and T^ex = 40 K are 
730 and 580 AU, respectively. Such numbers are comparable 
to the observed radius of the dust disk in scattered hght, about 
4.3" or 600 AU (Padgett et al., priv. comm., Schneider et 
al.,priv. comm.) 

To better understand the observed fine profile of '^CO, a 
simple disk model has been used together with a ray-tracing 
program (Hogerheijde & van der Tak 2000) to compute the line 
profile. Such models are better suited to constrain parameters 
using method with the observed spectrum. The populations 
of the ^^CO rotational levels are assumed to be in LTE through- 
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with the best model fit with a disk of i=2Q°, /f=600 AU. 
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Inl a nl jIn J|i,flhi 



5 _^ 10 
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15 



Fig. 5. '-CO 3-2, "CO 3-2 and C'^0 2-1 Unes toward GQ Lup, 
shifted by +4, +2 and K respectively. 



out the entire disk. The surface density of the disk is defined as 
S(r) = 2o(r/l AU)-^-^ with a typical value of 4.57 x lO^"^ cm'^ 
at 1 AU. The disk is assumed to be isothermal in the vertical di- 
rection. The abundance of CO is set to be 5x10"'' with respect 
to H2, except in areas with temperatures below 30 K where CO 
is frozen out onto the grains. There, abundances are taken to be 
5x10"'". It is assumed that the gas is in Keplerian orbits around 
a 0.5 Mo star. The inner radius is set at 0. 1 AU. The radial tem- 
perature profile is only constrained by a radial exponent of -0.5 
and the luminosity. The turbulent velocity is taken to be 0.1 km 
s-\ following Thi et al. (2001), Qi et al. (2004). 

The parameters to be fitted are the inclination, the outer ra- 
dius of the disk and the temperature at the inner radius. The 
best fit is presented in Fig 4. It has a best fit outer radius be- 
tween 400 and 700 AU, consistent with the result of 600 AU 
found above. The line profile is not very sensitive to the outer 
radius or the CO abundance in the outer areas because tem- 
peratures are too low to have a significant impact on the 3-2 
emission. The observed line is quite sensitive to the inclination 
of the disk, which is constrained to 20° + 5°. The temperature 
at the inner radius is taken to be 1200 K, but model results are 
not sensitive to the precise choice of temperature. 

4.2. GQ Lup 

The classical T Tauri star GQ Lup has recently attracted atten- 
tion because adaptive optics imaging has revealed a substel- 
lar companion, perhaps a massive Jupiter-like planet, with a 
separation of -100 AU (Neuhauser et al. 2005). As for all ob- 
jects in our sample, the GQ Lup spectrum shows an infrared 
excess and a silicate emission feature at 10 /^m, indicating the 
presence of a dust disk. This object is potentially interesting 
to further our understanding of the influence of a companion 
on the disk structure and evolution. High signal to noise ob- 
servations were taken toward this object (Fig. 5), but Uke most 
other sources in our sample only difliise intersteUar cloud ma- 



terial is detected in '^CO, as evidenced in the observations at 
ofi^-positions. The '^CO spectrum shows a hint of two peaks 
centered around yLSR=5 km s"' of equal height at rmb=0.25 
K, but deeper integrations are needed to confirm this. C'^O 2-1 
was not detected at all down to 5 mK (2cr) in a 0.25 km s"' bin. 
This constrains the maximum column density of C'^O directly 
associated with the source to ~ 2 x 10^^ cm"^. An unresolved 
disk would then be limited to a gas mass of < 2.3 x 10"^ M© 
for Ci'^O/Hz = 1.67 X 10"^ or < 2.3 x 10"^ Mq for a lower ratio 
of Ci^O/H2 = 1.67 X 10"^. The mass derived from the 1.3 mm 
flux of 38 mJy is 2.1x10""', assuming a gas/dust ratio of 100, 
^dust = 30 K and a dust opacity k of 0.01 cm^g"'(Nurnberger 
et al. 1997). 

5. Concluding remarks 

The main conclusions of this work are as foUows: 

- A '^CO /=3-2 survey of 21 classical T Tauri stars with 
dust disks in Lupus taken in frequency-switched mode has 
revealed mostly extended diffuse cloud material with a few 
magnitudes of extinction and complex velocity structure. 

- A C'**0 7=3-2 and 2-1 survey of the same 21 sources taken 
in beam-switched mode shows dense compact circumstel- 
lar material associated with only 2 sources, HT Lup and Sz 
73. 

- Only a single source, IM Lup (Sz 82), has been found with 
clear evidence for a gas-rich disk. Modeling of the Une in- 
tensity and line profile indicates a large disk with radius 
~600 AU and inclination 20°. This source will be a prime 
target for future interferometer observations in the southern 
sky, especially with the Atacama Large Millimeter Array. 

The lack of detection of gas disks in Lupus in the current sur- 
vey does not imply that the gas has been dissipated from these 
disks; rather, it shows the limitations of single-dish CO line 
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observations to reveal them. Future searches need to either fo- 
cus on high density single-dish tracers such as HCO"^ 7=4-3 
or HCN 7=4-3 to distinguish the warm, dense gas in the disk 
from the cold, diffuse cloud material or use newly commisioned 
interferometers such as the SubMillimeter Array (SMA). From 
comparison with other disks in Taurus (Thi et al. 2004, Greaves 
2004), these lines are expected to have peak intensities of at 
most 0.1 K in a single-dish 15" beam and will be very difficult 
to detect. Increased sensitivity to compact emission is only pos- 
sible with interferometer studies. Interferometric surveys will 
therefore be essential to determine the fraction of gas-rich disks 
in Lupus and thus the timescale for gas dissipation in this star 
forming region compared with other regions such as Taurus. 
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